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Abstract. We study the internal dynamics of the rich galaxy cluster ABGC 209 on 
the basis of new spectroscopic and photometric data. The distribution in redshift 
shows that ABCG 209 is a well isolated peak of 112 detected member galaxies at 
z — 0.209, characterised by a high value of the line-of-sight velocity dispersion, 
(Jv = 1250-1400 km s^^, on the whole observed area (1 h~^ Mpc from the cluster 
center) , that leads to a virial mass of M = 1.6-2.2 X lO^'^h"^ M0 within the virial 
radius, assuming the dynamical equilibrium. The presence of a velocity gradient in 
the velocity field, the elongation in the spatial distribution of the colour-selected 
likely cluster members, the elongation of the X-ray contour levels in the Chandra 
image, and the elongation of cD galaxy show that ABCG 209 is characterised by a 
preferential NW-SE direction. We also find a significant deviation of the velocity 
distribution from a Gaussian, and relevant evidence of substructure and dynamical 
segregation. All these facts show that ABCG 209 is a strongly evolving cluster, 
possibly in an advanced phase of merging. 

Key words. Galaxies: clusters: general - Galaxies: clusters: individual: ABCG 209 
- Galaxies: distances and redshifts - intergalactic medium - Cosmology: observa- 
tions 
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1. Introduction 



The investigation of clusters of galaxies offers a rare possibility to link many aspects of 
astrophysics and cosmology and, in particular, to understand the processes that lead to 
the formation of structures. 

In hierarchical clustering cosmological scenarios galaxy clusters form from the accre- 
tions of subunits. Numerical simulations show that clusters form preferentially through 



anisotropic accretion of subclusters along filaments (West et al. 1991; Katz & White 1993; 



Cen & Ostriker 1994; Colberg et al 1998, 199!;). The signature of this anisotropic cluster 
formation is the cluster elongation along the main accretion filaments (e.g., Roettiger 



et al. 1997). Therefore the knowledge of the properties of galaxy clusters, plays an im- 
portant role in the study of large-scale structure (LSS) formation and in constraining 
cosmological models. 

By studying the structure of galaxy clusters it is possible to discriminate between 



different cosmological models (e.g., Richstone et al. 1992; Lacey & Cole 1994; Thomas et 



al. 199q ). In fact, in a low-density universe the clustering tends to freeze at (z+1) = fl^ , 
while in a high-density universe it continues to grow to the present day. This implies that 
clusters in a low-density universe are expected to be dynamically more relaxed and to 
have less subsystems, called substructures. 

On small scales, clusters appear as complex systems involving a variety of interacting 
components (galaxies. X-ray emitting gas, dark matter). A large fraction of clusters 
(30%-40%) contain substructures, as shown by optical and X-ray studies (e.g., Baier & 



Ziener |1977| ; Geller & Beers |1982| ; Girardi et al. |1997| ; Jones & Forman |1999|) and by recent 
results coming from the gravitational lensing effect (e.g., Athreya et al. 2002| ; Dahle et 



al. 2002), suggesting that they are still in the dynamical relaxation phase. Indeed, there 
is growing evidence that these subsystems arise from merging of groups and/or clusters 
(cf. Buote 2002; and Girardi & Biviano 2002| for reviews). 

Very recently, it was also suggested that the presence of radio halos and relics in 
clusters is indicative of a cluster merger. Merger shocks, with velocities larger than 10'^ 
km s~^, convert a fraction of the shock energy into acceleration of pre-existing relativistic 
particles and provide the large amount of energy necessary for magnetic field amplification 
(Feretti 2000[ ). This mechanism has been proposed to explain the radio halos and relics 



in clusters (Brunetti et al. 2001). 

The properties of the brightest cluster members (BCMs) are related to the cluster 
merger. Most BCMs are located very close to the center of the parent cluster. In many 
cases the major axis of the BCM is aligned along the major axis of the cluster and of 
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the surrounding LSS (e.g., Binggeli [l98^ ; Dantas et al. |1997| ; Durret et al. |1998[ ). These 
properties can be explained if BCMs form by coalescence of the central brightest galaxies 



of the merging subclusters (Johnstone et al. 1991). 

The optical spectroscopy of member galaxies is the most powerful tool to investigate 
the dynamics of cluster mergers, since it provides direct information on the cluster ve- 
locity field. However this is often an ardue investigation due to the limited number of 
galaxies usually available to trace the internal cluster velocity. To date, at medium and 
high redshifts (z ^ 0.2), only few clusters are really well sampled in the velocity space 



(with > 100 members; e.g., Carlberg et al 1996; Czoske et al. 2002). 

In order to gain insight into the physics of the cluster formation, we carried out a 
spectroscopic and photometric study of the cluster ABCG209, at Z'-~^0.2 (Kristian et al. 



197S; Wilkinson & Oke 1978; Fetisova 1981), which is a rich, very X-ray luminous and 



hot cluster (richness class R = 3, Abell et al. |1989[ Lx(0.1 - 2.4 keV) ~ 14 K,^ 10^^^ erg 
s"\ Ebehng et al. |1996| ; Tx 10 keV, Rizza et al. |1998| ) . The first evidence for its 
complex dynamical status came from the significant irregularity in the X-ray emission 
with a trimodal peak (Rizza et al. 1998| ). Moreover, Giovannini et al. (1999) have recently 
found evidence for the presence of a possible extended Radio-emission. 

The paper is organised as follows. In Sect. 2 we present the new spectroscopic data 
and the data reduction. The derivation of the redshifts is presented in Sect. 3. In Sect. 4 
we analyse the dynamics of the cluster, and in Sect. 5 we complete the dynamical analysis 
with the information coming either from optical imaging or from X-ray data. In Sect. 6 
we discuss the results in terms of two pictures of the dynamical status of ABCG209. 
Finally, a summary of the main results is given in Sect. 7. 

Throughout the paper, we assume a flat cosmology with JIm = 0.3 and f^A = 0.7. For 
the sake of simplicity in rescaling, we adopt a Hubble constant of 100 h km s^^ Mpc^^. 
In this assumption, 1 arcmin corresponds to ^ 0.144 Mpc. Unless otherwise stated, we 
give errors at the 68% confidence level (hereafter c.l.). 



2. Observations and data reduction 

The data were collected at the ESQ New Technology Telescope (NTT) with the ESQ 
Multi Mode Instrument (EMMI) in October 2001. 

2.1. Spectroscopy 

Spectroscopic data have been obtained with the multi-object spectroscopy (MOS) 
mode of EMMI. Targets were randomly selected by using preliminary R-band images 
(Toxp=180 s) to construct the multislit plates. We acquired five masks in four fields 
(field of view 5' x 8.6'), with different position angles on the sky, allocating a total of 166 
slitlets (alligment stars included). In order to better sample the denser cluster region, we 
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covered this region with two masks and with the overlap of the other three masks. We 
exposed the masks with integration times from 0.75 to 3 hr with the EMMI-Grism#2, 
yielding a dispersion of 2.8 A/pix and a resolution of 8 AFWHM, in the spectral 
range 385 - 900 nm. 

Each scientific exposure (as well as flat fields and calibration lamps) was bias sub- 
tracted. The individual spectra were extracted and flat field corrected. Cosmic rays were 
rejected in two steps. First, we removed the cosmic rays lying close to the objects by 
interpolation between adjacent pixels, then we combined the different exposures by us- 
ing the IRAjj] task IMCOMBINE with the algorithm CRREJECT (the positions of the 
objects in different exposures were checked before). Wavelength calibration was obtained 
using He-Ar lamp spectra. The typical r.m.s. scatter around the dispersion relation was 
~ 15 km s""'^ . The positions of the objects in the slits were defined interactively using 
the IRAF package APEXTRACT. The exact object position within the slit was traced in 
the dispersion direction and fitted with a low order polynomial to allow for atmospheric 
refraction. The spectra were then sky subtracted and the rows containing the object were 
averaged to produce the one-dimensional spectra. 

The signal-to-noise ratio per pixel of the one-dimensional galaxy spectra ranges from 
about 5 to 20 in the region 380-500 nm. 



2.2. Photometry 

A field of 9.2' x 8.6' (1.2 x I.l h^^ Mpc^ at z=0.209) was observed in B-, V- and R- 
bands pointed to the center of the cluster. Two additional adjacent fields were observed in 
V-band in order to sample the cluster at large distance from the center (out to ~ 1.5 h^^ 
Mpc). The reduction of the photometric data will be described elsewhere (Mercurio et 



al. 2002, in preparation). In order to derive magnitudes and colours we used the software 



SExtractor (Bertin & Arnouts 1996) to measure the Kron magnitude (Kron 198C) in 
an adaptive aperture equal to a ■ rx, where rK is the Kron radius and a is a constant. 



Following Bertin & Arnouts (1996), we chose a = 2.5, for which it is expected that the 
Kron magnitude encloses ~ 94% of the total fiux of the source. We use the photometric 
data to investigate the colour segregation in Sect. 5. 



3. Redshifts measurements 



Redshifts were derived using the cross-correlation technique (Tonry & Davis 1981 ), as im 



plemented in the RVSAO package. We adopted galaxy spectral templates from Kennicutt 



^ IRAF is distributed by the National Optical Astronomy Observatories, which are operated 
by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 



A. Mercurio et al.: Internal Dynamics of ABCG209 at z~0.21 5 

( |1992D , corresponding to morphological types EL, SO, Sa, Sb, Sc, Ir. The correlation was 
computed in the Fourier domain. 

Out of the 166 spectra, 112 turned out to be cluster members (seven of which observed 
twice), 22 are stars, 8 are nearby galaxies, 1 is foreground and 6 are background galaxies. 
In 10 cases we could not determine the redshift. 

In order to estimate the uncertainties in the redshift measurements, we considered 
the error calculated with the cross-correlation technique, which is based on the width of 
the peak and on the amplitude of the antisymmetric noise in the cross correlation. The 
wavelength calibration errors (~ 15 km s^^ ) turned out to be negligible in this respect. 
The errors derived from the cross-correlation are however known to be smaller than the 
true errors (e.g., Malumuth et al. |1992| ; BardeUi et al. [l994 Quintana et al. ^000[ ). We 
checked the error estimates by comparing the redshifts computed for the seven galaxies 
for which we had duplicate measurements. The two data sets agree with one-to-one 
relation, but a reasonable value of for the fit was obtained when the errors derived 
from the cross-correlation were multiplied by a correction factor (^ 1.75). A similar 
correction was obtained by Malumuth et al. (1992; 1.6), Bardelli et al. ( 1994| ; 1.87), and 
Quintana et al. ( ^000 ; 1.57). External errors, which are however not relevant in the study 
of internal dynamics, cannot be estimated since only two previous redshifts are available 
for ABCG209. 

The catalogue of the spectroscopic sample is presented in Table 0, which includes: 
identification number of each galaxy, ID (Col. 1); right ascension and declination (J2000), 
a and S (Col. 2 and 3); V magnitude (Col. 4); B-R colour (Col. 5); heliocentric corrected 
redshift z (Col. 6) with the uncertainty Az (Col. 7); 

Our spectroscopic sample is 60% complete for V<20 mag and drops steeply to 30% 
completeness at V<21 mag; these levels are reached both for the external and the inter- 
nal cluster regions. The spatial distribution of galaxies with measured redshifts does not 
show any obvious global bias. Only around the brightest cluster member (cD galaxy; c.f. 



La Barbera et al. 2002 in preparation) the spatial coverage is less complete because of ge- 
ometrical restrictions. This does not affect our dynamical analysis because we investigate 
structures at scales larger than ~ 0.1 h~^ Mpc . 



4. Dynamical analysis 

4.1. Member selection and global properties 

ABCG 209 appears as a well isolated peak in the redshift space. The analysis of the 
velocity distribution based on the one-dimensional adaptive kernel technique (Pisani 



1993, as implemented by Fadda et al. 1996 and Girardi et al. 1996) confirms the existence 



of a single peak at z ~ 0.209. 
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Fig. 1. Distribution of redshifts of the cluster members. The big and small arrows indicate 
the mean cluster redshift and the redshift of the cD galaxy, respectively. 



Fig. |l| shows the redshift distribution of the 112 cluster members. The mean redshift 
in the present sample is < z >= 0.2090 ± 0.0004, as derived by the biweight estimator 



(Beers et al. |I990[) . 

In order to determine the cluster center, we applied the two-dimensional adaptive 
kernel technique to galaxy positions. The center of the most dense peak {a= 01 31 52.70, 
(5= -13 36 41.9) is close to the position of cD galaxy (a= 01 31 52.54, -13 36 40.4). 

We estimated the line-of-sight (LOS) velocity dispersion, a^, by using the biweight 
estimator (ROSTAT package; Beers et al. 1990| ). By applying the relativistic correction 



and the usual correction for velocity errors (Danese et al. 198C ) , we obtained CTv — 1394lgg 
km s^"'^, where errors were estimated with the bootstrap method. 

In order to check for possible variation of < z > and with increasing radius we plot 
the integrated mean velocity and velocity dispersion profiles in Fig. ^. The measure of 
mean redshift and velocity dispersion sharply vary in the internal cluster region although 
the large associated errors do not allow to claim for a statistically significant variation. 
On the other hand, in the external cluster regions, where the number of galaxies is larger, 
the estimates of < z > and Uv are quite robust. 

Assuming that the system is in dynamical equilibrium, the value of leads to a value 
of the radius of the collapsed, quasi-virialized region Rvir ~ 1.78 h^^ Mpc (cf. Eq. (1) 
of Girardi & Mezzetti 2001). Therefore our spectroscopic data sample about half of the 
virial region {'^ 1 h^^ Mpc). Under the same assumption, we estimated the mass of the 
system by applying the virial method. In particular, for the surface term correction to the 
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Fig. 2. Integrated mean redshift and LOS velocity dispersion profiles (upper and lower 
panel, respectively), where the mean and dispersion at a given (projected) radius from 
the cluster center is estimated by considering all galaxies within that radius. The error 
bars at the 68% c.l. are shown. In the lower panel, the faint line gives the profile after the 
rejection of seven possible interlopers according to the "shifting gapper" method (see text 
in Sect. 4.2), and the horizontal lines represent X-ray temperature with the respective 
errors (cf. Sect. 5.3) transformed in CTv imposing /3spoc = 1, where /3spec — fv/(kT/^mp), 
with fjL the mean molecular weight and nip the proton mass. 



standard virial mass we assumed a value of 20% (cf. Girardi & Mezzetti 2001), obtaining 
M(< Rvir) = 2.25l°;^^ x IQi^ h-^ Mq . 



4.2. Possible contamination effects 



We further explored the reliability of CTv related to the possibility of contamination by 
interlopers. 
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Fig. 3. Velocity in the cluster rest frame vs. (projected) clustercentric distance for the 112 
selected members. The standard application of the "shifting gapper" method would reject 
galaxies as cluster members indicated by squares, and a more conservative application 
also galaxies indicated by crosses. 



The cluster rest-frame velocity vs. projected clustercentric distance is shown in Fig. 
Although no obvious case of outliers is present, we applied the procedure of the "shifting 



gapper" by Fadda et al. (199£). This procedure rejects galaxies that are too far in velocity 
from the main body of galaxies of the whole cluster within a fixed bin, shifting along the 
distance from the cluster center. According to the prescriptions in Fadda et al. ( 1996| ), 
we used a gap of 1000 km s^^ and a bin of 0.4 h^^ Mpc, or large enough to include at 
least 15 galaxies. 

In this way we rejected three galaxies (indicated by squares in Fig. However, for 
ABCG 209 the results are too much sensitive to small changes of the adopted parameters. 
For instance, with a bin of 0.5 h~^ Mpc no galaxy was rejected, while seven galaxies 
were rejected with a bin of 0.3 h^^ Mpc (cf. crosses in Fig. ^). In the last case, we 
obtained a value of ~ 1250l;gg km s^^, which is slightly smaller than- (although fully 
consistent with-) the value computed in Sect. 4.1 (cf. also the velocity dispersion profile 
in Fig. ^). With this value of the velocity dispersion, the computation for the mass 
within Rvir = 1.59 h^^ Mpc leads to a value of total mass M(< Rvii) = l-62io;4g x lO^"^ 
h-i M0 . 

We determined the galaxy density and the integrated LOS velocity dispersion along 
the sequence of galaxies with decreasing density beginning with the cluster center (Kittler 



1976; Pisani 1996). As shown in Fig. H the density profile has only two minor peaks, 



possibly due to non complete sampling, and their rejection leads to a small variation in 



A. Mercurio et al.: Internal Dynamics of ABCG209 at z~0.21 



9 



1.5 



1 - 



0.5 - 








J L 



J L 



20 40 60 80 
R / [h-> Mpc] 



100 



Fig. 4. Integrated LOS velocity dispersion profile CTv (solid line, in units of 10^ km s"'^) 
and galaxy number density p (dashed line, in arbitrary units) along the sequences be- 
ginning with the center of the cluster. The thinner line gives the (Tv-profile after the 
rejection of galaxies belonging to the two small density peaks. 



the estimate of CTv and in the behaviour of velocity dispersion profile (cf. with Figure 2 



of Girardi et al. 1996 where the large effect caused by a close system is shown) . 

We conclude that the contamination by interlopers cannot explain the high value 
of velocity dispersion, which is therefore connected to the peculiarity of the internal 
dynamics of the cluster itself. In fact, a high value of CTv is already present in the internal 
cluster region, namely within ~ 0.2-0.3 h^^ Mpc (cf. Fig. ||), where the contamination 
is expected to be negligible. 



4.3. Velocity distribution 

In order to better characterise the velocity distribution, we considered its kurtosis K , 
skewness S, scaled tail index ST I, and the probability associated to the W-test P{W), 



(cf. Shapiro & Wilk |1965D . 

For the kurtosis and the skewness we found K — 2.61 ± 0.45, S — —0.20 ± 0.23, 
respectively, i.e. values that are consistent with a Gaussian distribution (reference value 
K—i^ 5=0). The ST I indicates the amount of the elongation in a sample relative to 
the standard Gaussian. This is an alternative to the classical kurtosis estimator, based 
on the data set percentiles as determined from the order statistics (e.g., Rosenberger & 



Gasko 1983). The Gaussian, which is by definition neutrally elongated, has ST I = 1.0. 



Heavier-tailed distributions have STI ^ 1.25 (see, e.g.. Beers et al. 1991). For our data 
STI = 1.254, indicating a heavier-tailed distribution, with 95%-99% c.l. (cf. Table 2 of 
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Fig. 5. Distribution of member galaxies separated into three clumps according the one- 
dimensional KMM test. The plot is centered on the cluster center. Dots, open circles, 
and crosses indicate low-, intermediate-, and high-velocity clumps, respectively. 



Bird & Beers 1993 ). On the other hand, the W-test rejected the null hypothesis of a 
Gaussian parent distribution, with only a marginal significance at 92.3% c.l. 

In order to detect possible subclumps in the velocity distribution, we applied the 



KMM algorithm (cf. Ashman et al 1994 and refs. therein). By taking the face value of 
maximum likelihood statistics, we found a marginal evidence that a mixture of three 
Gaussians (of ni = 13, 112 = 70, and ns = 29 members at mean redshift zi = 0.1988, 
Z2 = 0.2078, and Z3 = 0.2154) is a better description of velocity distribution (at 91.2% 
C.L). For the clumps we estimated a velocity dispersion of a^i = 337 km s^^ , (Tv2 = 668 
km s^^ , and (7^3 — 545 km s^^ . Fig. ^ shows the spatial distributions of the three 
clumps. The second and the third clumps are clearly spatially segregated. 



4.4. Velocity field 

The above result suggested to investigate the velocity field in more detail. To this aim, 
we divided galaxies into a low- and a high-velocity sample (with respect to the mean 
redshift). Fig. || shows that low- and high- velocity galaxies are clearly segregated in 
the E-W direction, the two distributions being different at the 99.88% c.l., according 
to the two-dimensional Kolmogorov-Smirnov test (hereafter 2DKS-test; cf. Fasano & 



Franceschini 1987, as implemented by Press et al. 1992). 

We then looked for a possible velocity gradient by means of a multi-linear fit (e.g.. 



implemented by NAG Fortran Workstation Handbook. 1986) to the observed velocity 
field (cf. also den Hartog & Katgert 1996| ; Girardi et al. 1996). We found marginal ev- 
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Fig. 6. Spatial distribution of member galaxies: the larger the circle, the smaller is the 
radial velocity. The figure is centered on the cluster center. The solid and the thin lines 
indicate the position angle of the cluster velocity gradient and relative errors, respectively. 



idence (c.l. 95.2%) for the presence of a velocity gradient in the direction SE-NW, at 
position angle PA= Wllgy degrees (cf. Fig. In order to assess the significance of the 
velocity gradient, we performed 1000 Monte Carlo simulations by randomly shuffling the 
galaxy velocities and determined the coefficient of multiple determination {R"^) for each 
of them. We then defined the significance of the velocity gradient as the fraction of times 
in which the of the simulation was smaller than the observed R^ . 



4.5. 3D substructure analysis 

In order to check for the presence of substructure, we combined velocity and position 



information by computing the A~statistics p devised by Dressier & Schectman (1988) . 
We found a value of 162 for the A parameter, which gives the cumulative deviation of 
the local kinematical parameters (velocity mean and velocity dispersion) from the global 
cluster parameters. The significance of substructure was checked by running 1000 Monte 
Carlo simulations, randomly shuffling the galaxy velocities, obtaining a significance level 
of 98.7%. 

This indicates that the cluster has a complex structure. 

The technique by Dressier & Schectman does not allow a direct identification of 
galaxies belonging to the detected substructure; however it can roughly identify the 



^ For each galaxy, the deviation 5 is defined as 5^ = (ll/(7^)[(vi — v) + {a\ ~ a) ], where 
subscript 1 denotes the average over the 10 neighbours of the galaxy. A is the sum of the 5 of 
the individual galaxies. 
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Fig. 7. Spatial distribution on the sky of the 112 cluster members, each marked by a 
circle: the larger the circle, the larger is the deviation 8 of the local parameters from the 
global cluster parameters. The figure shows evidence for substructure according to the 
Dressier & Schectman test. The plot is centered on the cluster center. 

positions of substructures. To this aim, in Fig. ^the galaxies are marked by circles whose 
diameter is proportional to the deviation (5 of the individual parameters (position and 
velocity) from the mean cluster parameters. 

A group of galaxies with high velocity is the likely cause of large values of 5 in the 
external East cluster region. The other possible substructure lies in the well sampled 
cluster region, closer than 1 arcmin SW to the cluster center . 



5. Galaxy populations and the hot gas 

5.1. Luminosity and colour segregation 

In order to unravel a possible luminosity segregation, we divided the sample in a low 
and a high-luminosity subsamples by using the median V-magnitude= 19.46 (53 and 54 
galaxies, respectively), and applied them the standard means test and F-test (Press et 



al. 1992). We found no significant difference between the two samples. 

We also looked for possible colour segregation, by dividing the sample in a blue and a 
red subsamples relative to the median colour B-R= 2.32 of the spectroscopic sample (41 
and 44 galaxies, respectively). There is a slight differences in the peaks of the velocity 
distributions, Zbiuc = 0.2076 and z^^d = 0.2095, which gives a marginal probability of 
difference according to the means test (93.5%). The velocity dispersions are different 
being (Tv.biuo = 1462j:i45 km s^^ larger than (Tv.rcd = 99318^ km s^^ (at the 98.6% c.L, 
according to the F-test), see also Fig. Moreover, the two subsamples differ in the 
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Fig. 8. Integrated LOS velocity dispersion profile as in Fig. |^. The solid and dashed lines 
give the profiles of redder and bluer galaxies, respectively. The faint lines give the profiles 
of bright and faint galaxies (upper and lower lines, respectively). 
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Fig. 9. Spatial distribution of the 85 member galaxies having B-R colour information: 
circles and crosses give redder and bluer galaxies, respectively. The plot is centered on 
the cluster center. 



distribution of galaxy positions (97.4% according to the 2DKS test), with blue galaxies 
lying mainly in the SW~region (cf. Fig. ^). 

The reddest galaxies (B-R> 2.43) are characterised by a still smaller velocity dis- 
persion (732 km s~^) and the bluest (B-R< 2.19) by a larger velocity dispersion (1689 
km s^^). Fig. |l^ plots velocities vs. clustercentric distance for the four quartiles of the 
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Fig. 10. Velocity in the cluster rest-frame vs. (projected) clustercentric distance for the 
85 members having B-R colour information. Large circles, small circles, small crosses, and 
large crosses indicate galaxies belonging to the four quartiles of the colour distribution, 
from reddest to bluest galaxies. 



colour distribution: the reddest galaxies appear well concentrated around the mean clus- 
ter velocity. From a more quantitative point of view, there is a significant negative corre- 
lation between the B~R colour and the absolute value of the velocity in the cluster rest 
frame jvif | at the 99.3% c.l. (for a Kendall correlation coefficient of —0.18). 

The velocity of the cD galaxy {z— 0.2095) shows no evidence of peculiarity according 



to the Indicator test by Gebhardt & Beers (1991). Moreover, the cD galaxy shows an 
elongation in the NE-SW direction (cf. Fig. ^ see La Barbera et al. 2002, in preparation, 
for more details) 



5.2. Substructures with photometric data 

We performed a two-dimensional analysis to detect subclumps using the photometric 
sample of galaxies for which we have B-R colours. This allows us to take advantage of 
the larger size of the photometric sample compared to the spectroscopic one. Galaxies 
were selected within ±0.5 mag of the B-R vs. R colour-magnitude relation and within 
the completeness limit magnitude R= 22 mag. The colour-magnitude relation was de- 
termined on the N=85 spectroscopically confirmed cluster members, for which we have 
B-R colours. The above selection leads to 3.68 < (B-R) -|- 0.097 R < 4.68, with a sample 
size of N=392 galaxies. 

Fig. |l^ shows the B-R vs. R colour-magnitude diagram for the 596 galaxies detected 
in the R-band in the ABCG 209 central field. The sample includes also the 85 member 
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Fig. 11. V-band image of cD galaxy (0.7' x 0.5'), with North at top and East to left. 



galaxies with known redshift. The colour-magnitude relation was derived by using a 
linear regression based on the biweight estimator (see Beers et al. 1990). 

Fig. |l^ - upper panel - shows the projected galaxy distribution is clustered and 
elongated in the SE-NW direction. Notice that there is no clear clump of galaxies per- 
fectly centered on the cD galaxy, whose position is coincident with the center determined 
in Sect. 4.1 from redshift data only. In order to investigate this question, we noted that 
galaxies with redshift data are brighter than 22 R-mag, being the ^ 90% of them brighter 
than 19.5 mag. We give the distributions of 110 galaxies and 282 galaxies brighter and 
fainter than 19.5 mag, respectively, in central and lower panels (Fig. Brighter galax- 
ies are centered around the cD, while fainter galaxies show some clumps aligned in the 
SE-NW direction; in particular the main clump. Eastern with respect the cD galaxies, 
coincides with the secondary peak found in our analysis of Chandra X-ray data (see 
below) . 



5.3. Analysis of X-ray data 

The X-ray analysis of ABCG 209 was performed by using the Chandra ACIS-I obser- 
vation 800030 (exposure ID#522) stored in the Chandra archive^. The pointing has an 
exposure time of 10 ksec. 

Data reduction was performed by using the package CIAO^Chandra Interactive 
Analysis of Observations). First, we removed events from the level 2 event list with 



http:/ /asc. harvard.edu/cda/ 



* CIAO is freely available at ittp://asc. harvard.edu/ciao/ 
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Fig. 12. B-R vs. R colour-magnitude diagram for galaxies in the central field of 
ABCG209. Filled circles: spectroscopically confirmed cluster members; open circles: all 
galaxies within the completeness limit. 



a status not equal to zero and with grades one, five and seven. Then, we selected all 
events with energy between 0.3 and 10 keV. In addition, we cleaned bad offsets and ex- 
amined the data on a chip by chip basis, filtering out bad columns and removing times 
when the count rate exceeded three standard deviations from the mean count rate per 
3.3 second interval. We then cleaned each chip for flickering pixels, i.e. times where a 
pixel had events in two sequential 3.3 second intervals. The resulting exposure time for 
the reduced data is 9.86 ks. 

The temperature of the ICM was computed extracting the spectrum of the cluster 
within a circular aperture of 3 arcminutes radius around the cluster center. Freezing 
the absorbing galactic hydrogen column density at 1.64 10^° cm~^, computed from the 



HI maps by Dickey & Lockman (1990), we fitted a Raymond-Smith spectrum using 
the CIAO package Sherpa with a statistics. We found a best fitting temperature of 
Tx = 10.21^;^ keV. 

In order to detect possible substructures in ABCG 209 we ran the task 
CIAO/Wavdetect on a subimage of 1600 by 1600 pixels containing the cluster. The task 
was run on different scales in order to search for substructure with different sizes. The 
significance threshold was set at 10^^. 
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arcmin 





Fig. 13. Spatial distribution on the sky and relative isodensity contour map of likely 
cluster members selected on the basis of the B-R colour-magnitude relation, recovered 
by adaptive-kernel method (cf. Pisani |1993| , |1996D . The plots are centered on the cluster 
center. Top panel: the 392 selected galaxies (see text). The 1 arcmin circle is centered on 
the cD galaxy. Central panel: bright (R< 19.5) galaxies. The 1 arcmin circle is again cen- 
tered on the cD galaxy. Bottom panel: fainter (R> 19.5) galaxies. The 1 arcmin circle is 
centered on the secondary X-ray peak (cf. Sect. 5.3). In the figures the isodensity contours 
are snaced bv 1 • 10~^ galaxies /arcsec^ in the ton Danel and by 5 • 10~^ galaxies /arcsec^ 
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Fig. 14. Optical image of ABCG 209 with, superimposed, the X-ray contour levels of 
the Chandra image, with North at top and East to left. The ellipses represent the two 
substructures detected by Wavdetect. The primary X-ray peak coincides with the po- 
sition of the cD galaxy and thus with the peak of the bright galaxy distribution (cf. 
Fig. |l^ middle panel). The secondary X-ray peak coincides with the main peak of the 
faint galaxy distribution (cf. Fig. |l^ lower panel). 

The results are shown in Fig. Two ellipses are plotted representing two significant 
substructures found by Wavdetect. The principal one, located at a= 01 31 52.7 and S= 
-13 36 41, is centered on the cD galaxy, the left one is a secondary structure located at 
a= 01 31 55.7 and 6= -13 36 54, about 50 arcseconds (-120 h^^ kpc) East of the cD (cf. 
also Fig. [l^) . The secondary clump is well coincident with the Eastern clump detected by 



Rizza et al. (1998) by using ROSAT HRI X-ray data, while we do not find any significant 



substructure corresponding to their Western excess. 



6. Formation and evolution of ABCG 209 



The value we obtained for the LOS velocity dispersion ~ 1400 km s ^ (Sect. 4.1) is 
high compared to the typical, although similar values are sometimes found in clusters 



at intermediate redshifts (cf. Fadda et al. 1996; Mazure et al. 1996; Girardi & Mezzetti 



2001). In Sect. 4.2 we show that this high value of is not due to obvious interlopers: 
also a very restrictive application of the "shifting gapper" or the rejection of minor peaks 
in galaxy density lead to only a slightly smaller value of CTv — 1250 km s~^. The high 
global value of CTv ~ 1250-1400 km s~^ is consistent with the high value of Tx ~ 10 
keV coming from the X-ray analysis (which, assuming /3spcc — 1, would correspond to 
CTv ~ 1300 km s-i , cf. Fig. |) and with the high value of Lx(0.1 - 2.4 keV) - 14 h^^ lO"*" 



erg s ^ (Ebeling et al. |l996| ; cf. with Lx-ctv relation by e.g., Wu et al. |l999[ Girardi & 
Mezzetti ^001]) . 

Therefore, on the basis of the global properties only, one could assume that ABCG 209 
is not far from dynamical equilibrium and rely on large virial mass estimate here 
computed M(< 1.78 h^^ Mpc) = 2.25^0.11 x 10^^ M(< 1.59 h^^ Mpc) = 

l.62to\tl X 10^5 h-i Mq . 

On the other hand, the analysis of the integral velocity dispersion profile shows that a 
high value of is already reached within the central cluster region of 0.2-0.3 h^^ Mpc. 
This suggests the possibility that a mix of clumps at different mean velocity causes the 
high value of the velocity dispersion. 

A deeper analysis shows that ABCG 209 is currently undergoing a dynamical evolu- 
tion. We find evidence for a preferential SE-NW direction as indicated by: a) the presence 
of a velocity gradient; b) the elongation in the spatial distribution of the colour-selected 



A. Mercurio et al.: Internal Dynamics ol ABCG209 at z~0.21 19 

likely cluster members; c) the elongation of the X-ray contour levels in the Chandra 
image; d) the elongation of the cD galaxy. 

In particular, velocity gradients are rarely found in clusters (e.g., den Hartog & 
Katgert 1996D and could be produced by rotation, by presence of internal substructures, 
and by presence of other structures on larger scales such as nearby clusters, surrounding 



superclusters, or filaments (e.g.. West 1994; Praton & Schneider 1994). 

The elongation of the cD galaxies, aligned along the major axis of the cluster and of 



the surrounding LSS (e.g., Binggeli 1982; Dantas et al. 1997; Durret et al. 1998), can be 
explained if BCMs form by coalescence of the central brightest galaxies of the merging 



subclusters (Johnstone et al. 1991). 

Other evidence that this cluster is far from dynamical equilibrium comes out from 
deviation of velocity distribution from Gaussian, spatial and kincmatical segregation 
of members with different B-R colour, and evidence of substructure as given by the 



Dressler-Schectman test ( 1988 ). In particular, although a difference in Cv between blue 
and red members is common in all clusters (e.g., Carlberg et al. 1997), a displacement in 
mean velocity or in position center is more probably associated with a situation of non 



equilibrium (e.g., Bruzendorf & Meusinger 1999). 

Possible subclumps, if any, cannot be easily separated by using 112 cluster members 
only. The separation in three clumps by using one-dimensional KMM test and the visual 
inspection of Dressler-Schectman map of Fig. ^ represent only two tentative approaches 
(cf. Sects. 4.3 and 4.5). 

A further step in the detection of possible subclumps is achieved by using the two- 
dimensional complete sample of all ~ 400 colour-selected cluster members. In fact, several 
clumps appear well aligned in the SE-NW direction, and, in particular, the main clump 
revealed by faint galaxies (R> 19.5 mag) coincides with the secondary X-ray peak as 



found in our analysis of Chandra data (cf. also Rizza et al. 1998 ). This result suggests 
that one or several minor clumps are merging along the SE-NW direction, with a main 
clump hosting the cD galaxy. In particular, the presence of a significant velocity gradient 
suggests that the merging plane is not parallel to the plane of the sky. 

The presence of a secondary clump as a well distinct unit in X-ray data suggests that 
this clump and maybe the whole cluster might be in a pre-merger phase. 

The strong luminosity segregation found for colour-selected galaxies could however 
suggests a possible alternative dynamical scenario. In fact, very appealingly, galaxies of 
different luminosity could trace the dynamics of cluster mergers in a different way. A first 



evidence was given by Biviano et al. (1996): they found that the two central dominant 
galaxies of the Coma cluster are surrounded by luminous galaxies, accompanied by the 
two main X-ray peaks, while the distribution of faint galaxies tend to form a structure not 
centered with one of the two dominant galaxies, but rather coincident with a secondary 
peak detected in X-ray. The observational scenario of ABCG 209 has some similarities 
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with the situation of Coma: bright galaxies (R< 19.5 mag) are concentrated around the 
cD galaxy, which coincides with the primary X-ray peak, while faint galaxies (R> 19.5 
mag) show several peaks, the main of which. Eastern with respect the cD galaxy, is 
confirmed by the position of the secondary X-ray peak. Therefore, following Biviano et 
al., we might speculate that the merging is in an advanced phase, where faint galaxies 
trace the forming structure of the cluster, while more luminous galaxies still trace the 
remnant of the core-halo structure of a pre-merging clump, which could be so dense 
to survive for a long time after the merging (as suggested by numerical simulations 



Gonzalez-Casado et al. 1994). An extended Radio-emission would support an advanced 
merging phase, but its presence is still uncertain due to the existence of strong discrete 



sources (Giovannini et al. 1999). The comparison of present Radio-image with our galaxy 
distribution shows that the diffuse source is located around the central cD galaxy with 
an extension toward North. 

Unfortunately, redshift data are available only for luminous galaxies and we cannot 
investigate the nature of subclumps inferred from the two-dimensional distribution of 
faint galaxies. At the same time, available X-ray data are not deep enough to look for 
possible variations of temperature in the region of interest. 



7. Summary 

In order to study the internal dynamics of the rich galaxy cluster ABCG 209, we obtained 
spectra for 159 objects in the cluster region based on MOS observations carried out at 
the ESO New Technology Telescope. Out of these spectra, we analysed 119 galaxies: 112 
turn out to be cluster members, 1 is foreground and 6 are background galaxies. 

ABCG 209 appears as a well isolated peak in the redshift distribution centered at 
z — 0.209, characterized by a very high value of the LOS velocity dispersion: = 1250- 
1400 km s"\ that results in a virial mass of M = 1.6-2.2 x lO^^h^^ Mq within Rvir. 
The analysis of the velocity dispersion profile show that such high value of CTv is already 
reached in the central cluster region (< 0.2-0.3 h^^ Mpc). 

The main results of the present study may be summarised as follows. 

— ABCG 209 is characterised by a preferential SE-NW direction as indicated by: a) the 
presence of a velocity gradient in the velocity field; b) the elongation in the spatial 
distribution of colour-selected likely cluster members; c) the elongation of the X-ray 
contour levels in the Chandra image; d) the elongation of the cD galaxy. 

— We find significant deviations of velocity distribution from Gaussian. 

— Red and blue members are spatially and kinematically segregated. 

— There is significant evidence of substructure, as shown by the Dressier & Schectman 
test. 
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— The two dimensional distribution of the colour-selected likely members shows a 
strong luminosity segregation: bright galaxies R < 19.5 are centered around the cD 
galaxy, while faint galaxies R > 19.5 show some clumps. The main one, Eastern with 
respect to the cD galaxy, is well coincident with the secondary X-ray peak. 

This observational scenario suggests that ABCG 209 is presently undergoing strong 
dynamical evolution. Present results suggest the merging of two or more subclumps 
along the SE-NW direction in a plane which is not parallel to the plane of sky, but 
cannot discriminate between two alternative pictures. The merging might be in a very 
early dynamical status, where clumps are still in the pre-merging phase, or in a more 
advanced status, where luminous galaxies trace the remnant of the core-halo structure 
of a pre-merging clump hosting the cD galaxy. 

The connection between the dynamics and the properties of galaxy populations will 
be discussed in a forthcoming paper. 
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Table 1. Spectroscopic data. 



ID a 5 Y B-R z Az 



1 


01 


31 


33.81 


-13 


32 


22.9 


18.32 




0.2191 


0.0003 


2 


01 


31 


33.82 


-13 


38 


28.5 


19.27 


2.43 


0.2075 


0.0002 


3 


01 


31 


33.86 


-13 


32 


43.4 


19.46 




0.1998 


0.0003 


4 


01 


31 


34.11 


-13 


32 


26.8 


19.90 




0.2000 


0.0004 


5 


01 


31 


34.26 


-13 


38 


13.1 


20.17 


2.24 


0.2145 


0.0005 


6 


01 


31 


35.37 


-13 


31 


18.9 


17.36 




0.2068 


0.0004 


7 


01 


31 


35.37 


-13 


37 


17.4 


18.58 


2.38 


0.2087 


0.0002 


8 


01 


31 


35.42 


-13 


34 


52.0 


19.10 


2.37 


0.2073 


0.0002 


9 


01 


31 


35.61 


-13 


32 


51.3 


20.17 




0.2090 


0.0002 


10 


01 


31 


36.51 


-13 


33 


45.2 


19.83 


2.36 


0.2072 


0.0004 


11* 


01 


31 


36.54 


-13 


37 


56.1 


19.89 


1.37 


0.2620 


0.0005 


12 


01 


31 


36.76 


-13 


33 


26.1 


21.10 


2.21 


0.2064 


0.0004 


13 


01 


31 


36.87 


-13 


39 


19.3 


19.73 


2.26 


0.2039 


0.0002 


14a 


01 


31 


37.19 


-13 


30 


04.3 


20.65 




0.3650 


0.0003 


15 


01 


31 


37.23 


-13 


32 


09.9 


20.03 




0.2084 


0.0003 


16 


01 


31 


37.38 


-13 


34 


46.3 


18.71 


2.35 


0.2134 


0.0003 


IT" 


01 


31 


37.99 


-13 


36 


01.6 


21.38 


2.67 


0.3997 


0.0004 


18 


01 


31 


38.06 


-13 


33 


35.9 


19.87 


2.49 


0.2073 


0.0002 


19 


01 


31 


39.06 


-13 


33 


39.9 


19.64 


2.33 


0.2120 


0.0002 


20 


01 


31 


38.69 


-13 


35 


54.7 


20.11 


2.17 


0.2084 


0.0003 


21 


01 


31 


39.89 


-13 


35 


45.1 


18.85 


1.68 


0.2073 


0.0002 


22 


01 


31 


40.05 


-13 


32 


08.9 


19.51 




0.2096 


0.0003 


23 


01 


31 


40.17 


-13 


36 


06.2 


19.18 


2.18 


0.1995 


0.0002 


24 


01 


31 


40.58 


-13 


36 


32.9 


19.30 


2.18 


0.2052 


0.0003 


25 


01 


31 


40.76 


-13 


34 


17.0 


18.58 


2.33 


0.2051 


0.0002 


26 


01 


31 


40.94 


-13 


37 


36.5 


19.10 


1.66 


0.2181 


0.0008 


27 
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Table 1. Continued. 
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Table 1. Continued. 



ID 




a 




5 




V 


B-R 


z 




87 


01 


oz 


Ui.oZ 


-13 


33 


34.0 


18.60 


2.52 


0.2083 


0.0003 


88 


01 


oz 




-13 


36 


15.7 




2.33 


0.2131 


0.0004 


89 


01 


oz 


ui.yi 


-13 


35 


31.2 


18.65 




0.1984 


0.0002 


90 


01 


oz 


no 1 s 

Uz.lo 


-13 


35 


51.0 


19.06 


2.10 


0.2091 


0.0003 


91 


01 


Q9 
oZ 


no Al 


-13 


32 


13.0 


19.90 




0.1982 


0.0004 


92 


01 


oZ 


no 0/1 


-13 


39 


39.4 


19.48 


2.48 


0.2110 


0.0002 


93 


01 


oz 


n'^ n7 


-13 


36 


44.0 


19.95 


2.33 


0.2102 


0.0005 


94 


01 


oz 


uo.yo 


-13 


35 


54.4 


19.92 


2.33 


0.2137 


0.0003 


95 


01 


Q9 
oZ 


nQ 07 


-13 


38 


01.2 


20.31 


2.16 


0.2116 


0.0003 


96 


01 


Q9 
oZ 


nQ c;n 


-13 


31 


59.1 


19.65 




0.1971 


0.0003 


97 


01 


Q9 
OZ 


C\A OQ 


-13 


39 


53.9 


17.82 


2.55 


0.2125 


0.0002 


98 


01 


Q9 
OZ 


n/1 Qc; 


-13 


37 


26.3 


18.36 


2.50 


0.2061 


0.0002 


99 


01 


Q9 
OZ 


nn 


-13 


40 


35.1 


19.77 


2.34 


0.2132 


0.0002 


100 


01 


Q9 
OZ 


nci no 


-13 


33 


42.0 


20.34 


2.30 


0.2087 


0.0003 


101 


01 


Q9 
OZ 


nc; 


-13 


37 


33.3 


19.28 


2.35 


0.2074 


0.0002 


102'' 


01 


Q9 
oZ 


n7 OQ 

u ( -zy 


-13 


37 


30.3 


19.06 


1.63 


0.1744 


0.0006 


103 


01 


Q9 
oZ 


n7 OS 


-13 


41 


31.4 


19.86 


2.23 


0.1997 


0.0002 


104 


01 


oZ 


1 n 11 


-13 


37 


24.3 


19.54 




0.2161 


0.0002 


105 


01 


oZ 


1 Q/1 


-13 


34 


21.1 


19.74 




0.2134 


0.0004 


106 


01 


oz 


1 f^l 
IZ.Oi 


-13 


41 


18.1 






0.2145 


0.0003 


107 


01 


oZ 


1 1 


-13 


39 


46.0 






0.2087 


0.0003 


108^ 


01 


^9 
oZ 


1 ^ 00 

lo.zz 


-13 


31 


09.1 


22.34 




0.2597 


0.0005 


109 


01 


Q9 
OZ 




-13 


38 


08.5 


17.90 




0.2175 


0.0002 


110 


01 


Q9 
OZ 


1 /I c;o 


-13 


32 


29.1 


19.39 




0.2145 


0.0002 


111 


01 


Q9 
OZ 


1 nn 


-13 


41 


13.9 






0.2158 


0.0003 


112 


01 


Q9 
OZ 


iO.OD 


-13 


37 


49.2 


18.94 




0.2149 


0.0002 


113 


01 


32 


15.84 


-13 


35 


41.0 


20.71 




0.2124 


0.0004 


114 


01 


32 


16.00 


-13 


38 


34.6 


19.57 




0.2203 


0.0003 


115 


01 


32 


16.20 


-13 


32 


38.8 


18.97 




0.2034 


0.0003 


116 


01 


32 


16.98 


-13 


33 


01.1 


20.64 




0.2117 


0.0004 


117 


01 


32 


16.76 


-13 


35 


03.7 


20.81 




0.2138 


0.0004 


118 


01 


32 


17.26 


-13 


38 


06.1 


19.74 




0.2105 


0.0003 


119 


01 


32 


17.48 


-13 


38 


42.3 


19.84 




0.1963 


0.0006 



Background galaxy. 
^ Foreground galaxy. 



This figure "isofote.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0209536v3 



